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The running charm-quark mass in the MS scheme is determined from weighted finite energy 
QCD sum rules (FESR) involving the vector current correlator. Only the short distance expansion 
of this correlator is used, together with integration kernels (weights) involving positive powers of s, 
the squared energy. The optimal kernels are found to be a simple pinched kernel, and polynomials 
of the Legendre type. The former kernel reduces potential duality violations near the real axis 
in the complex s-plane, and the latter allows to extend the analysis to energy regions beyond the 
end point of the data. These kernels, together with the high energy expansion of the correlator, 
weigh the experimental and theoretical information differently from e.g. inverse moments FESR. 
Current, state of the art results for the vector correlator up to four-loop order in perturbative QCD 
are used in the FESR, together with the latest experimental data. The integration in the complex 
s-plane is performed using three different methods, fixed order perturbation theory (FOPT), 
contour improved perturbation theory (CIPT), and a fixed renormalization scale p (FMUPT). The 
final result is m c (3GeV) = 1008 ± 26MeV, in a wide region of stability against changes in the 
integration radius so in the complex s-plane. 
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Considerable progress has been made over the years to 
extract accurate values of the quark masses by con- 
fronting QCD with experimental data in the framework 
of QCD sum rules We consider in this paper the 
case of the charm-quark mass. Following pioneering ap- 
proaches , it has been customary to write the Operator 
Product Expansion (OPE) of the vector correlator with 
the scale invariant mass m c (m c ) as the expansion param- 
eter. This is in contrast to the case of light quarks where 
the square of the four-momentum Q 2 = —q 2 is the nat- 
ural large expansion parameter. If one were to consider 
the charm quark mass in e.g. the MS scheme at a typical 
scale of 3 GeV, then the fact that m c (3 GeV) ~ 1 GeV 
should be a matter of some concern regarding this ex- 
pansion [3]. In any case, the latest determinations based 
on inverse moment QCD sum rules claim an accuracy at 
the 1% level These inverse moments require QCD 
knowledge of the vector correlator in the low energy re- 
gion, around the open charm threshold, as well as in the 
high energy region. An alternative approach involving 
positive moment sum rules, requiring QCD information 
only at high energy, and naturally suited to determine 
fhc(n) in the MS scheme, was proposed some time ago in 
[j| . We follow this approach here, and make use of state 



of the art QCD information up to four-loop level 
together with the latest experimental data [U 
consider Finite Energy Sum Rules (FESR) weighted by 
two types of integration kernels. The first type is the 
so called pinched kernel, which has been shown to min- 
imize potential duality violations close to the real axis 
in the complex s-plane [22[. The second type is a poly- 
nomial (Legendre) kernel tuned to minimize the impact 
of the energy region where the data is either poor or 
non-existent. Effectively, these Legendre-type integra- 
tion kernels, involving several positive powers of the en- 
ergy, allow for an extension of the analysis beyond the 
end point of the data. Both types of kernels have been 
used successfully in the light quark sector to study the 
saturation of chiral sum rules [23j . to quantify duality 
violations '24S\ , to extract the values of the vacuum con- 
densates in the OPE (2||, and to determine the quark 
condensates [26| and light quark masses [27j . 
We begin by considering the vector current correlator 
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where V fl (x) = c(x)7^c(x). Invoking Cauchy's theorem 
in the complex s-plane (— q 2 = Q 2 = s) one has 

/•so i i r 
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where p(s) is an arbitrary but analytical integration ker- 
nel, and 



Imll(s) = -L R c (s) 
iZn 
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with R c (s) the standard R-ratio for charm production. 
The perturbative QCD (PQCD) expression of II(s) at 
high energies can be written as 

n( s )|p Q ci^£(^YV)( s ), (4) 
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where 
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and to c stands for the running charm-quark mass in the 
A/5-scheme. The complete analytical PQCD result for 
H(s)\pqcd up to order O [a 2 {m 2 /s) & ] is given in @, and 
exact results for IIq 3 '' and LI^ are from Q. The func- 
tion Tl^ is known exactly up to a constant Q which 
has been estimated using Pade approximants in At 
five-loop order 0(a 4 ) the full logarithmic terms for li^ 
are given in [loj . and for LI^ in 11 1. Since there is in- 
complete knowledge at this loop-order, we shall use the 
available information as a measure of the truncation er- 
ror in PQCD. 

The fundamental parameters entering the QCD corre- 
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FIG. 1: The running mass m c (3 GeV) in the MS scheme for a 
fixed fi = 3.0 GeV (FMUPT), as a function of so and using the 
pinched kernel Eq. (6). Results from FOPT and CIPT lead 
to similar stability regions. Errors are due to uncertainties 
in the data, to the values of a s and the gluon condensate, to 
PQCD truncation, and to changes in fi of ±50%. 

lator are the running strong coupling a s (/i 2 ), and the 
gluon condensate. Regarding the strong coupling, we 



use the latest comprehensive update analysis at the 
r-scale [H[ which gives a s (M 2 ) = 0.342 ± 0.012, or 
a s (M 2 ) = 0.1213 ± 0.0014. This value agrees within 
errors with a recent determination [29j | from e + e~ anni- 
hilation data: a„{M%) = 0.1172 ± 0.0051. However, a 
world average with much smaller errors [30l ] gives the re- 
sult a s (M|) = 0.1184± 0.0007, in agreement with lattice 
QCD results [3lj . We shall include both sets of values of 
a s in our determination of m c . In the non-perturbative 
sector the leading power correction in the OPE involves 
the gluon condensate, i.e. ((a s /n)G 2 ). The latest value 
of the gluon condensate [25[ , extracted from the ALEPH 
data on r-decays, is {{a 8 /ir)G 2 ) = (0.046± 0.012) GeV 4 , 
for A QCD = 300 MeV, and <(a s /7r)G 2 ) = (0.021 ± 
0.006) GeV 4 for A QCD = 350 MeV, where A QCD stands 
for the QCD scale in the MS'-scheme. While the gluon 
condensate is renormalization group invariant, there is an 
unavoidable dependence on a s when its value is extracted 
from data using QCD sum rules. In fact, the sum rules in- 
volve the difference between hadronic data integrals and 
PQCD integrals, with the latter being strongly depen- 
dent on the value of a s . In other words, the uncertainty 
in a s induces an uncertainty in the condensate. Ex- 
trapolating the above results to include a value Aqcd = 
380 MeV, in line with the latest determination of a s (M 2 ) 
[Hj], leads to {(a s /n)G 2 ) = (0.01 ± 0.01) GeV 4 . This 
large uncertainty in the value of the gluon condensate 
will have a very small impact on our results for fh c . This 
is because the high energy expansion of the vector corre- 
lator, together with integration kernels involving positive 
powers of s, tend to reduce considerably the importance 
of this term. This is in contrast to the inverse moments 
method which enhances this contribution with increasing 
inverse powers of s. Something similar happens with the 
impact of the uncertainty in a s . Results for m c from in- 
verse moments are far more sensitive to a s than results 
from sum rules involving the high energy behaviour of 
the correlator and positive powers of s, as with pinched 
or Legendre integration kernels. 

Turning to the experimental data, our analysis follows 
closely that of 0], [32j |. For the first two resonances 
we use the latest data from the Particle Data Group 
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3.096916(11) GeV, T JH , = 5.55(14) keV, 



3.68609(4) GeV, T 



M2s) = 2.35(4) keV. The 
first two resonances are followed by the open charm re- 
gion where the contribution from the light quark sector 
R u ds must be subtracted as background from the to- 
tal R-ratio Rtot- This we do as in [32j. In the region 
3.97 GeV < y/s < 4.26 GeV we only use CLEO data 
[2l| as they have the least error. Regarding the three 
data sets from BES [li|-[2(|, we assume conservatively 
that the systematic uncertainties are not entirely inde- 
pendent and add them linearly, rather than in quadra- 
ture, but we treat them as independent from the CLEO 
data set [21(, and thus add these in quadrature. The re- 
gion s = 25 — 49 GeV 2 has no data, beyond which there 



is CLEO data up to s c 
compatible with PQCD. 



110 GeV z . The latter is fully 
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We discuss next the integration kernels p(s) in Eq.(2), 
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FIG. 2: The running mass m c (3 GeV) in the MS scheme as 
a function of so, using the Legendre-type polynomial V5 with 
si = 24 GeV 2 . Polynomials of different order lead to equally 
wide stability regions. The errors are due to uncertainties 
in the data, the values of a s and the gluon condensate, and 
changes in fi of ± 50%. 

and first introduce the pinched kernel 

p(s) = 1 - — . (6) 
so 

As shown in the past [22j - [HJ, this kernel suppresses 
potential duality violations close to the real axis in the 
complex s-plane. Clearly, there is a large variety of al- 
ternative functional forms vanishing at s = sq. Having 
considered many of these alternatives, the simple choice 
above turns out to be the optimal. In any case, higher 
powers of s will be part of the Legendre-type kernels we 
discuss next. The purpose of these kernels is to extend 
the analysis beyond the end point of the data, as achieved 
e.g. in the analysis of duality violations using the ALEPH 
r-decay data |2J] . In the present case, there is no data in 
the interval 25 GeV 2 < s < 50 GeV 2 , while the data for 
s > 50 GeV 2 agrees with PQCD. Hence, we introduce a 
kernel p(s) such that 

p(s)=V n [x(s)}, (7) 

where 

= 2s-(s + s 1 ) 
s - si 

with so > 81, and V n (x) are the standard Legendre poly- 
nomials, i.e. Vi(x) = x, V%{x) — (5a; 3 — 3x)/2, etc., 
which satisfy the constraint 

f ° s k T n [x{s)} ds = , (9) 

J si 

where si ~ 24 GeV 2 , and so varies in the region where 
there is no data. In the hadronic sector, the Legendre- 
type kernels provide extra weight to the well known reso- 
nance region on account of their rapid growth for s < s±. 
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FOPT 


996 
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11 


25 


FMUPT 


988 


9 
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25 


CIPT 


983 


9 


1 


2 


1 


16 


25 



TARLE I: Results for m c (3 GeV) (in McV) in the A4S scheme using 
the pinched kernel Eq.(6) in FOPT, CIPT, and in FMUPT using a fixed 
scale fi — 3 GeV. Listed results are at sq — 17 GeV 2 , a representative 
value inside the wide stability region (see Fig. 1). The uncertainties arc 
due to the data (Exp.), to a s (Aa 8 ), to truncation of PQCD (Trunc), 
to the gluon condensate (NP), and to variation inside the stability re- 
gion in so (Var.). All uncertainties arc added in quadrature except for 
that in so which is conservatively added linearly. The FMUPT total 
error includes an uncertainty in of ± 50%. 
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v 3 


100 


1008 


24 
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25 




160 


1008 


24 
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26 


r 5 


200 


1007 


22 


1 


8 


2 


23 




300 


1008 


23 


1 


9 


2 


25 



TABLE II: Results for m c (3 GeV) (in McV) in the MS scheme ob- 
tained using various Legendre polynomial-type kernels, Eqs. (8)-(10), 
with si — 24 GeV 2 , fj, = 3 GeV, and representative values of so in the 
remarkably wide stability region (see Fig. 2). The uncertainties are 
due to the data (Exp.), to a s (Ack s ), to changes in fi by ±50% (A/^), 
and to the gluon condensate (NP), which arc added in quadrature. The 
variation in the stability region is negligible, and so is the truncation 
uncertainty. 



In Fig. 1 we show the result for to c (3 GeV) as a function 
of so using a fixed scale fj, = 3 GeV (FMUPT). The errors 
shown are the result of adding in quadrature the uncer- 
tainties due to experiment, to the values of a s and the 
gluon condensate, and to the truncation of PQCD which 
is taken as the difference between the known four-loop 
result and the partially known five-loop expression (to 
zeroth order in m c ). This FMUPT total error includes 
an uncertainty in /j, of ± 50% (not present in FOPT or 
CIPT). This combined error is then conservatively added 
linearly to the uncertainty due to variation of so in the 
stability region. Results from FOPT and CIPT are simi- 
lar, and fully compatible within errors. Numerical values 
are given in Table I. In Fig. 2 we show m c (3GeV) as 
a function of so for the kernel V5 with si = 24 GeV 2 . 
Legendre-type polynomials of other order lead to basi- 
cally the same results, which are very stable in a remark- 
ably wide region. Numerical values are listed in Table 
II where the uncertainties are due to the data (Exp.), 
to a s (Aa s ), to changes in // by ±50% (A/x), and to 
the gluon condensate (NP). Variation within the stabil- 
ity region gives a negligible uncertainty, and so does the 
PQCD truncation error. In comparison with values in 
Table I, one notices the larger impact of experimental 
uncertainties. This is because so is much larger, thus in- 
cluding more data in the hadronic integral. Since quark 
mass terms in PQCD are suppressed by inverse powers 
of s, the Legendre-type kernels are expected to improve 
convergence for increasing energies. We choose as our 
final result the one from the Legendre-type kernels due 
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to its stability against changes in sq, and use the pinched 
kernel values in Table I as confirmation. In this case 

m c (3 GcV) = 1008 ± 26 MeV , (10) 

in good agreement within errors with the result from in- 
verse moments QCD sum rules [4|, other recent deter- 
minations 32], [34} . as well as lattice QCD [3l|. Trans- 
lated into a scale invariant mass, the above result gives 
w c (tti c ) = 1319 ± 26 MeV for the value used here for 
the strong coupling, a s (M|) = 0.1213 ± 0.0014. Using 
instead a s (M|) = 0.1189 ± 0.0020, as in 0, changes this 
mass to fh c {m c ) = 1299 ± 26 MeV. 
In closing we briefly discuss two convergence issues, (a) 
the convergence pattern in a s of the perturbative QCD 
integral as a function of m Cl and (b) the convergence pat- 
tern of the values of m c resulting from the FESR trun- 
cated at different orders in a s . In the first case the per- 
turbative contour integral 

z( So ) = -J-i" P ( s )n( s )ds (11) 
2ni Mm) 

is a function of both the mass ifi and the cou- 
pling a s . One can investigate the convergence of 
this integral as a function of a s for different val- 
ues of ra c . Using a typical integration kernel, e.g. 
p(s) = 7^4 and the representative value m = 1 GeV, 



we find = 1.4176 GeV 2 , 1^ = 1.4563 GeV 2 , 

J( 2 ) = 1.4541 GeV 2 , 1^ = 1.4512 GeV 2 and 

= 1.4524 GeV 2 , where the upper index in 
indicates the power of a s . If one were to use precisely 
the same input parameters, but with a higher quark 
mass, e.g. m = 5 GeV, the convergence would expect- 
edly be lost. To restore it one only needs to increase 
correspondingly the scale fi. Turning to the second 
issue, and using again the typical kernel 'P4, one can 
calculate m c (3 GeV) truncated at different orders in a s 
as a further check on the convergence. The results are 
fh [ c ] = 986 MeV, rh^ } = 1011 MeV, m^ 2) = 1010 MeV, 
mc ] = 1008 MeV, and m£ 4) = 1009 MeV, where the 
upper index in mP" indicates the power of a s . Consid- 
ering the difference between fh^ and fh^ would give 
an uncertainty of roughly 1 MeV, which is very much 
smaller than that from the 50% variation in the scale fi 
considered in our analysis. 
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